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We report the results of first-principles electronic-structure calculations for su- 
perconducting BaeCeo- Unlike the A^Cqq superconductors, this new compound 
shows strong Ba-C hybridization in the valence and conduction regions, mixed co- 
valent/ionic bonding character, partial charge transfer, and insulating zero-gap band 
structure. 
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Most of the known superconducting Ceo fullerides exhibit the stoichiometry A^Cqo or 
A 2 BCqq (A, B= alkali), form fee lattices, and are known (or believed) to exhibit essentially 
complete charge transfer from alkali to Ceo- Broadly speaking, all of these compounds 
exhibit similar electronic structure, in which the half-filled valence band (VB) is formed 
primarily from C orbitals of p z symmetry. Superconductivity has also been observed in 5:1 
solid solutions of Ca in Ceo [0- In contrast to the A 3 Cqo compounds, theoretical studies of 
hypothetical (stoichiometric) Ca n Ceo indicate substantial hybridization of C and Ca orbitals 
in the VB regime |§, but the relevance of these results to the real (non-stoichiometric) 
material is difficult to gauge. 

Recently, synthesis of a crystalline, stoichiometric, Ba-intercalated fulleride supercon- 
ductor, Ba 6 Ceo, has been reported, with T c =7 K ||. Naive electron counting, assuming 
divalent Ba and rigid K 6 C 60 -like bandstructure, would suggest this new compound to be 
a 0.2-eV gap insulator consisting of very highly charged C^ 12 molecules. Here, we present 
first-principles electronic structure results for BaeCgo that definitively predict strong Ba-C 
hybridization near the Fermi level, incomplete charge transfer and, surprisingly, a zero- gap 
insulating band structure. 

Like the alkali A 6 Cq compounds, Ba 6 C 60 forms a bee Bravais lattice, with two-fold 
planar crosses of Ba ions decorating the faces of the conventional unit cell. The aspect ratio 
of the crosses is identical to that in K 6 Cgo 0, and the lattice constant of 11.171 A is 2% 
smaller than for K 6 Cgo- As with K 6 Cgo, Rietveld refinement of x-ray powder patterns for 
Ba 6 C 60 indicate rotationally ordered C 60 molecules. 

We have studied the electronic structure and energetics of crystalline Ba 6 C 60 using the 
local-density approximation (LDA) to density-functional theory. The computational meth- 
ods were identical to those used in previous studies of K 6 C 60 and K 3 C 60 H^H, with a few 
modifications. In our local-orbital method 0, basis functions are expanded on a set of 
fourteen Gaussian exponents contracted into four s-type and three p-type functions for C, 
and 21 exponents contracted into seven s-type, five p-type, and four cf-type functions for 
Ba. The smallest Gaussian exponents were 0.24 for C and 0.15 for Ba; these basis sets give 
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eigenvalue spectra converged to a tolerance of ~0.02 eV for core, valence, and low-lying 
conduction bands. The charge density and potential are completely general and without 
shape approximation, and were iterated to self-consistency using the T point. All 696 elec- 
trons were explicitly included in the calculation and were treated on equal footing, without 
the need for pseudopotentials. The Ceperly- Alder exchange-correlation functional was used. 
Lowest-order relativistic effects were accounted for perturbatively, as described below. 

The resulting band structure and density of states (DOS) are shown in Fig. |l|, plotted 
along the same high-symmetry lines as for our previous study of KqCqq (cf. Fig. 2 of Ref. 
H). Although the correspondence between the two is clear, substantial differences arise from 
strong hybridization of C 2p and Ba 5d wave functions in the valence and conduction regions. 
Both the ti u band (filled in K 6 C 60 ) and the t lg band (empty in K 6 C 60 ) are filled in Ba 6 C 60 . 
Surprisingly, the electronic structure is again insulating, since there are no partially filled 
bands and the Fermi-level DOS is indeed zero. The magnitude of the gap is zero, to within 
the limits of accuracy of our calculation. This is strikingly at odds with the experimental 
finding of superconductivity (and therefore a metallic normal state) in BagCgo- Moreover, 
we note that while the LDA is known generally to give semiconductor band gaps smaller 
than experiment, we are unaware of any instance for which the LDA band structure of a 
known metal is insulating. We have exhaustively tested the quality of our basis set: for a 
variety of different sets of Gaussian exponents, the longest-ranged of which was 0.12 for C 
and 0.10 for Ba, the VB and CB varied by at most 0.02 eV, and by less than this near the 
Fermi level. 

A enlarged view of the DOS near the Fermi level is shown in Fig. ^|, along with atom- 
and symmetry-projected partial DOS. States within 1 eV of E-p have approximately 75% 
C 2p character and 25% Ba 5d character, with a very small admixture of C s character in 
the low-lying CB. The contribution from Ba s and p basis functions to bands in this energy 
region is negligible. 

We have also simulated relativistic corrections to the bands, by appending to the Hamil- 
tonian a non-local operator diagonal in the atomic-orbital Bloch basis. These atomic shifts 
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were taken over directly from differences between relativistic and non-relativistic atomic 
eigenvalues for Ba (the shifts for C are negligible). For Ba core states, the shifts are large 
and resulted in purely rigid shifts of the corresponding bands. For the 5d states, which 
contribute to the VB and lower CB, we used the degeneracy- weighted average of the atomic 
5d 3 /2 and 5<i 5 / 2 shifts ||, giving a single shift of +0.44 eV. Since the upper VB and lower CB 
states have comparable Ba 5d admixture, diagonalization of this new Hamiltonian resulted 
in nearly rigid band shifts in this energy region. Indeed, within 1 eV of Ep, the changes were 
at most 0.1 eV, and the zero-gap insulating character of the band structure was unchanged 
(the relativistically corrected bands are shown in two of the panels of Fig. p], and the partial 
DOS in Fig. ^includes this correction). Relativistic shifts for the Ba 6s atomic orbitals are 
quite small (—0.16 eV), and since these states contribute only to bands several eV above E F , 
these shifts have a negligible effect on the occupied eigenvalue spectrum and charge density. 

We have not systematically attempted to include spin-orbit (SO) effects here. The atomic 
SO splitting of the Ba 5d 3 / 2 and 5d 5 / 2 states is 0.11 eV 0. To a first approximation, one 
expects the resultant band splitting also to be proportional to the amount of 5d character, 
and so we estimate the SO splitting in both the upper VB and lower CB to be of order 0.05 
eV or less. Shifts of ±0.05 eV to E-p would result in a Fermi-level DOS, N(E-p), of order 
1-2 states/eV-cell-spin. This is 3-6 times lower than what we have previously calculated for 
K 3 C6o 0- If the pairing mechanism is the same in both materials and the coupling strengths 
are comparable, then such values for N(E-p) are far too low to account for the measured T c of 
7 K. This suggests a number of possible scenarios: (1) If SO splitting is indeed responsible for 
making the bands metallic, then either the BagCgo superconducting mechanism itself or the 
pairing strength (or both) is quite different from the standard fullerene superconductors. (2) 
If the one-electron (mean-field) approximation of the LDA does not apply to Ba 6 Cgo, then it 
is conceivable that a full accounting of electron-electron correlation leads to a metallic normal 
state. We view this as unlikely, since the standard result of treating electron correlation 
more accurately is to increase the band gap, or even to give an insulating ground state 
when the single-particle ground state is metallic. (3) In principle, it is possible that a self- 
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consistent spin-polarized calculation would lead to a metallic ground state. Such magnetic 
instabilities generally have as prerequisite a large value for N(Ep) in the paramagnetic (spin- 
unpolarized) state. Since we find a paramagnetic solution with N(E-p) ~ 0, we view this 
scenario as unlikely. (4) If we take seriously the zero-gap insulating state, then the possibility 
also exists that BCS-like pairing occurs not between states on the Fermi surface, but rather 
between states whose energies lie in some small interval E-p ± AE. Clearly, this picture 
would be favored by the existence of an extremely small, or zero, gap. Interestingly, this 
scenario also requires some minimum "threshold" coupling strength, below which pairing 
does not occur (in contrast to pairing of states on the Fermi surface, for which any non-zero 
coupling suffices). 

We consider now the question of charge transfer from Ba to Cgo- In a previous study 
of K 6 C 60 , we used Mulliken population analysis to compute the charge associated with 
each atom in the cell, and found essentially complete charge transfer from K to Ceo 0- 
For Ba 6 C 60 , Mulliken analysis associates 55.57 electrons with each Ba, suggesting a loss 
of only 0.43 electrons per atom. We caution, however, that Mulliken charges may not be 
meaningful when there is large wave-function overlap, and so we turn to direct methods. 
Integration of the total charge density within touching muffin-tin (MT) spheres, of radius 
R, gives the charges, Q(R), listed in the third column of Table |. Of course, these spheres 
do not represent the entire cell volume, nor are they the only plausible choice of integration 
volumes. Thus, the total integrated charges are less meaningful than the difference charges 
(relative to overlapping neutral atoms), AQ(R), shown in the fourth column. These values 
suggest that ~0.7 electrons are transferred from each Ba atom, in rough agreement with the 
Mulliken result. 

A more detailed description of the spatial distribution of Ba charge is given by the 
function AQ(r) for values of r less than the touching MT radii. For Ba, this function is 
shown in Fig. |^, for Rq < r < Rb&- It is clear that most of the 0.6 electron is lost from 
the region 1.6 A < r < 2.2 A. Moreover, we find that the Ba 5d and 6s atomic expectation 
values of r are 1.65 A and 2.33 A, respectively |§. These observations are consistent with a 
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picture in which charge is lost from the Ba 6s state, and donated partially to Ba 5c? states 
and partially to C p states. This description is also consistent with our prediction of the 
mixed hd/2p character of the highest-lying occupied bands. 

Some insight into the differences between the K- and Ba-intercalated fulleride crystals 
can be gained by examining the breakdown of an ionic model that was found to successfully 
describe K 6 Ceo- Since the details of this model have been described elsewhere, || we only 
present a brief qualitative discussion here. For a system of non-overlapping entitities, the 
charge transfer can be predicted by properly accounting for the Madelung contributions 
and the density-functional-based electron affinities (for the Ceo molecule) and the ionization 
energies (for the metallic dopants). When such a model is applied to K 6 C 60 , full charge 
transfer is expected at the experimental lattice constant of 11.39 A. Furthermore, since at 
this lattice constant the smallest C-K distance (3.2 A) is large compared to the sum of the 
ionic radii (2.1 A), corrections due to banding are expected to be small. Results from this 
ionic model were in excellent agreement with the charge transfer calculated self-consistently 
for K 6 C6o [§J- In contrast, when the same ionic model is applied to BagCgo, we find that 
the Ba atoms are expected to lose 1.2 electrons, roughly twice as much as calculated self- 
consistently in the crystal. This discrepancy can be understood by noting that a Ba atom 
in such a charge state would still have excess valence charge available for bonding with the 
C atoms. Further, for this charge state, the Ba 5d states exhibit a maximum at a radius 
of approximately 1.5-1.6 A and exhibit appreciable tails as far out as 2.6-2.8 A. Since 
the nearest-neighbor C-Ba distance is only 3.1 A, the partially filled Ba 5g? orbitals overlap 
strongly with the neighboring C 2p states, which allows for the formation of a covalent bond. 
This covalent bonding can be further strengthened by by allowing stronger overlap between 
the Ba and C states. At self-consistency, this is accomplished by reducing the ionicity of the 
Ba atoms over what is expected from the ionic model leading to longer ranged Ba d states. 
While the covalent bonds are strengthened, the Madelung stabilization decreases and the 
actual degree of charge transfer is arrived at by compromising between a purely covalent 
and purely ionic system. 
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To further characterize the degree of covalent vs. ionic bonding character, we compare 
the valence-electron density in BagCgo to that of KqCqq. In Fig. [|(a) we show the electron 
density from the filled t\ u band of K 6 C6o- The radially directed C-centered density lobes 
are plainly evident, and very little contribution from K states appears. (We also note a 
slight but definite polarization of charge toward the K + ions). In Fig. |](b) we show the 
density from the corresponding band of BaeCeo, and in Fig. f|(c) the t\ g density. A 
substantial amount of Ba-derived density is evident in each. Nearest-neighbor C-Ba pairs 
clearly show a substantial off-axis "bent" covalent bond, accompanied by strongly polarized 
back-bonding lobes. Weaker bonds are formed between next-nearest-neighbor C-Ba pairs. 
Given the substantial covalent character accompanied by (incomplete) charge transfer, we 
conclude that bonding in Ba 6 C6o is best described as mixed covalent /ionic. 

In summary, we have performed first-principles LDA calculations on superconducting 
BagCeo and find a number of unusual features. In contrast to the structurally identical 
compound KqCqq, for which we predicted full charge transfer and ionic bonding in, for 
Ba 6 Cgo we find only partial charge transfer and mixed covalent /ionic bonding. In contrast 
to the A 3 C 60 superconductors, we calculate a zero-gap insulating band structure, which 
appears to be at odds with the experimental finding of superconductivity. 

We thank W.E. Pickett and E.J. Mele for many illuminating discussions. This work was 
supported in part by the Laboratory for Research on the Structure of Matter (University of 
Pennsylvania) and in part by the Office of Naval Research. 
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FIGURES 

FIG. 1. Self-consistent electronic band structure for BaeCeo- The Fermi level is the energy zero. 
The dotted curves (for clarity shown only along the directions T — N — H) are the relativistically 
corrected bands described in the text. 

FIG. 2. Total DOS near the Fermi level (top panel); atom- and symmetry-projected partial 
DOS for carbon (middle panel) and barium (bottom panel). For the partial DOS, s components 
are shown as dotted curves and are magnified by a factor of 10 (100) for C (Ba); there are only 
negligible contributions from Ba p states in this energy region. Note also the factor of 3 scale 
difference between the C and Ba graphs. 

FIG. 3. Charge-difference function, AQ(r), in a sphere centered on a Ba atom. The arrows 
denote the expectation values of r for atomic Ba 5d and 6s states. 

FIG. 4. Valence-electron densities for (a) the t\ u band of K^Ceo, (b) the t\ u band of BagCeo, 
and (c) the t\ g band of BaeCgo- Each plotting plane contains four Cqo molecules (at the corners of 
the plot area) and four K + ions (marked by large crosses). The projected C positions are marked 
by small crosses, and the Cgo cage radius is shown by quarter circles. A heavy solid line marks the 
nearest-neighbor C-Ba internuclear axis. Adjacent contours are separated by 0.0005 a.u. 
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TABLES 

TABLE I. Integrated muffin-tin (MT) charges, Q(R), and charge differences, AQ(R), for the 
4 inequivalent atoms in the BaeCeo cell. C atoms are numbered by increasing distance from the 
nearest Ba atom. Charge differences are defined with respect to overlapping neutral atom electronic 
configurations. 



Atom MT radius, R (A) Q(R) AQ(R) 

Ci 0.70 4.05 +0.10 

C 2 0.70 4.04 +0.10 

C 2 0.70 4.06 +0.11 

Ba 2.37 59.24 -0.65 



10 



